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CYLINDRICAL  IMPLOSION  IN  SHALLOW  WATER  THEORY 
By  Cathleen  S.  Morawetz 

We  seek  here  a  theory  to  describe  the  behavior  of  a 
cylindrical  shock  in  shallow  water  as  it  contracts  to  the 
center  and  then  expands  again. 

The  corresponding  problem  for  gases  has  been  worked 
out  by  Guderley  [1]  for  strong  shocks.   Although  shallow 
water  theory  corresponds  to  a  gas  with  y  =  2,  the  cylindri- 
cal reflection  theory  breaks  down  because  the  shock  condi- 
tions are  essentially  different. 

The  solution  given  here  is  to  be  valid  only  in  the 
neighborhood  of  the  center  and  for  times  close  to  the  time 
when  the  shock  hits  the  center.   This  is  also  true  of  the 
reflection  solutions  in  gases.  However,  Guderley' s  solu- 
tion is  an  exact  solution  of  the  differential  equations 
corresponding  to  a  shock  of  prescribed  velocity,  which 
becomes  infinite  at  the  center,  impinging  on  a  region 
of  zero  pressure.   The  solution  here  for  a  shallow  water 
implosion  is  an  exact  solution  for  the  shallow  water 
approximation  only  if  the  water  ahead  of  the  shock  is  of 
zero  depth.   This  is  quite  unrealistic •   However,  for 
arbitrary  depth  ahead  the  given  solution  satisfies  the 
differential  equations  to  zero  order  in  (log  t)"  and 
the  boundary  conditions  corresponding  to  a  shock  velo- 
city which  becomes  infinite  In  a  prescribed  way,  at  the 
center.   It  is  expected  to  represent  the  flow  except  on 
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the  shock  paths . 

In  Section  1,  we  shall  find  the  differential  equations 
and  describe  the  method  of  finding  special  solutions.   In 
Section  2,  we  shall  determine  the  boundary  conditions  and 
find  the  solution  corresponding  to  zero  depth  ahead  of  the 
shock.   In  Section  3*  we  shall  modify  this  solution  and 
show  that  the  resultant  solution  satisfies  the  differen- 
tial equations  and  the  shock  conditions  near  the  center, 
up  to  terms  of  order  (log  t)~  .   In  Section  l\.t   we  discuss 
a  computed  flow, 

1 •   Special  Solutions 

According  to  shallow  water  theory  [2],  the  depth  Z 

and  the  horizontal  velocity  u  of  i^ater  satisfy  the 

equations, 

u.  +  u  u  +  gZ  =0 
t      r   e  r 


(1) 


Z,  +  u  Z„  +  Z  u  +  r^uZ  =  0 


if  the  flow  has  cylindrical  symmetry;   r  is  the  distance 
from  the  origin,   t  time,  and  g  the  acceleration  of  gra- 
vity. 

Setting  gZ  =  c  ,  equations  (1)  become 

u,  +  uu   +  2c  c   =0 
(2)  t     r      r 

^  -1 

2c,  +  2uc   +  cu  +  r  uc  =  0  . 

Note  that  c  has  the  dimensions  of  velocity  and  is 
the  speed  of  propagation  of  small  disturbances. 
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There  are  various  special  solutions  of  (2)  but  only  one 
class  need  be  considered  here.  All  others  lead  to  physically 
impossible  solutions.   We  reject  for  example  solutions  for 
which  Z  — >  0  as   t  — >  co   or  for  which  Z  — >  co   or 
u  — >  co   for   t  =  0,   r^  0  .  \le   must  permit  solutions  that 
are  singular  at  r  =  0,  t  =  0  . 

The  solutions  considered  are  the  so-called  progressive 
waves,  [3,  p«Lj.2lj.] 

(3)  u  =  5  U(,L) 


c  =  I  c(nj 


where 


Then     U     and     C      satisfy 

(5)  (1   -   U)rU.,   -    2Cr-,C„    +   U(U  -    1)    +   2C2   =   0 

w   'L  'l 

-   nUy,  +   2(1   -   U)C"1>lC,     +   I4.U  -    2   =   0 
or 

(6)  ((1   -   U)2   -   C2)^=  A 

((1  -   U)2   -  C2hG,{  =  BC 

where  A  =  U((l  -  U)2  -  2C2),   B  =  -  C2  -  (1  -  U)(|u  -  1)   or 
eliminating  Vt  we  get  a  nonlinear  ordinary  differential  equa- 
tion relating  U  and  C  : 

(7)  dC  =  BC  ' 

dU    A  ' 


I 


■ 


"" 
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Then  a  solution  of  (1)  is  given  by  (3)  if  C(U)   is  a 
solution  of  (7)  and  '1,  is  found  from  (6): 

(8)  to.=  ilz2^£  dD  . 

Prom  (8)  and  (L(.)  we  see,  as  in  the  case  of  gases  that 

2  2 

no  solution  may  cross   C   =  (1  -  U)    unless  A  =  0,  i.e. 

U  =  o,  C  =  +  ls  or  U  =  1,  C  =  0.   Otherwise  d*\,  changes 
sign  or  for  fixed  r,   dt  changes  sign  which  is  impossible. 

The  field  of  the  solutions  of  (7)  is  to  be  investiga- 
ted particularly  in  the  neighborhood  of  the  singular  points 
where  both  BC  and  A  vanish.   These  points  are 

U  =  1,   C  =  0 

u  =  i   c  =  +  -2- 

c  2/2 

U  =  0,   C  -  +  1 

U  =  0,   C  =  0  . 

First  from  (3),  to  have  a  solution  valid  on  the  line 
t  =  0  we  must  have  U  =  0,  C  =  0  for  If  I — >  0  .  Hence 
we  need  consider  only  solutions  of  (7)  which  pass  through 
(0,  0),   Prom  (7)  we  see  that  near   (0,  0) 

dC  • .  C 
"dU   U 

and  hence  U~  C  r^    constant  represent  solutions.   If  the 
solutions  are  to  be  continuous  across   t  =  0  then  the  slope 


. 


of  the  solution  curve  must  be  continuous  through  the  origin. 
Near  U  =  1,   C  =  0  we  find 

P   i 
dC  =  C(C^  +  2(1-U)) 

dU    2C2  -  (1-U)2 

By  the  theory  of  ordinary  nonlinear  differential  equations, 
this  singularity  is  a  node.   Except  for  the  solution  C  =  0, 
all  solutions  enter  U  =  1,   C  =  0   on  the  lines 
C  =  +  4(1-U),   To  determine  the  direction  of  increasing 

time,   di;>  0,  we  find  on  C2  =  i(l  -  U)   that  both 

P  P 

(1   -   U)      -  C        and     A      are  positive.      Hence      dlq  >  0     for 

HdU  >  0. 

The  point   (^  ,  +  -=— )   is  a  saddle  point  singularity, 
d  2/2 

The  singularity  at   (0,  +  1)   is  again  a  node. 

The  field  of  solutions  is  sketched  in  Figure  1,  cf. 

E3]  >   V*   ^4-26,  with  the  direction  of  increasing  time  indicated 

by  the  arrow, 
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2 •   Shock  Conditions  and  Corresponding  Special  Flow . 

The  conditions  satisfied  at  a  discontinuity  in  shallow 
x^ater  theory  are,  see  [2], 

(9)        Zq(uo  ~  s)  =  Z](ui  ~  ^)  (conservation  of  mass) 

Z  (u  -  S)   +  ~gZ  =  Z,  (u,-S)   +  ^gZ2  (conservation  of 
ov  o         2^o    lx  1       2to  j.  v 

momentum). 

Here  the  subscript   0  will  represent  the  state  ahead 
of  the  shock,  the  subscript   1   the  state  behind  the  shock. 
The  shock  velocity  is   S. 

A  third  condition  is  that  there  is  to  be  a  loss  of  energy 
across  the  shock.   This  is  equivalent  to  requiring 


u  -  S   >  C 

1  o     '     o 


or 


|u,  -  S|  <  C., 


9 

Setting  gZ  =  c   again,  we  have 


(10)  c2(uo  ~  S)  =  c2^  -  S) 


C2(uo  -  S)2  +  §c£  =c2(Ul  -  S)2  +|c^  . 


These  are  equivalent  to 


c 


2 


c2  =  -f   (  -  1  +  /L+  8(uo-S)2C;2) 

(ID  c2 

u,  -  S  =  -%   (u  -  S)  . 

cl 


■  •    ■  ,.  ■        ■ 


.' 


-7- 


The  suffixes  may  be  interchanged  to  determine  the  state 
behind  in  terms  of  the  state  ahead. 

In  the  case  u  =  0,  (11)  becomes 
"c2 


d2)  c2  =  -g  (-  i  +  J\  +  bs2c;2  ) 


Ux  =  S(l  -  c2c"2)  . 


In  the  limiting  case  where  c   — >  0  we  have 

1 
(13)  c^S  =  2^  c^2  S"l/2(1  +  0(co))— >  0 

a,/S  =  1  -  0(cQ)  — >  1   . 

We  note  that  we  are  led  to  the  same  result  if  we  fix 

c    and  let  S  — >  co  , 

o  ' 

(111.)  Cl/S   =  23Aj/2(l   +  0(S"1) 

u^/S   =  1   -   0(S"1)    . 

The  special  solution  which  corresponds  to  cavitation 
ahead,  c  =  0  ,  can  now  be  found  as  a  progressive  wave  in 
the  usual. way.   The  shock  is  given  by 


(1?)  dr  _  q  _  r 

dt  -  b  ~  t 

or 


r  =  kxt 
or 

-  v-1 
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The  shock  conditions  reduce  by  (13)>  (15>)  and  (3)   to 

(16)  U  =  1 

C  =  0 

We  take  any  solution  that  starts  from  U  =  1,  C  =  0 
in  the  direction  of  increasing  time  and  goes  through  U  =  0, 
C  =  0.   In  this  section  C  <  0,  since   t  <  0,  c  >  0  .   This 
solution  is  continued  then  through  the  origin.   It  is  to  be 
cut  off  by  a  reflected  shock  which  will  be  determined.   The 
state  behind  it  if  it  is  to  be  of  the  form  (3)rmust  be  the 
state  of  rest,   u  =  U  =  0  .   Otherwise  we  would  be  led  to 
dv]  <  0  or  cavitation. 

The  state  of  rest  is  given  by 

(17)  u  =  U  =  U2  =  0 


c  =  c2 


or  C  =  C0  =  CpT^ 

where   c~   is  to  be  determined. 

This  state  can  be  reached  from  states  of  the  form  (3) 
with  a  reflected  shock  given  by  r  =  kpt  provided  U  and 
C  satisfy,  by  (11)  and  (3), 

(18)  C2(k2)(U(k2)  -  1)  =  -  c|  k~2 

2  .  -2 
c~  k 


C2(k?)  =  -2^-2™  (-  1  +  /l  +  8c:2k2)  . 


2  ~2 
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These  values  of  U  and  C  lie  on  a  cur-®  passing  through 
(1,  0),  (0,  +  1),  and  extending  to   U  =  -  co  where  it 
approaches  the  lines  C  =  +  /2  asymptotically.   The  portion 
of  the  curve  where  |U  -  l|  >  G  ,  that  is  c?k~      >  1  corres- 
ponds to  a  state  ahead  of  a  shock. 

The  shock  curve  given  by  (18)  intersects  every  solution 
curve  passing  through  the  origin  into  U  <  0  .   To  every 
intersection  there  corresponds  a  value  of  c?k~   , 

We  are  still  left  with  a  large  family  of  solutions.   In 
fact  almost  every  solution  that  leaves  (1,  0)  and  enters 
(0,  0  -  )  is  possible  since  they  all  end  on  the  shock  curve. 
However,  the  solution  that  leads  to  the  lox^est  value  of 
Cpkp   --  that  is,  the  weakest  reflected  shock,  is  not  admiss- 
ible since  we  cannot  integrate  dvw .   This  solution  leaves 
(1,  0)  and  passes  through  the  saddle  point  {-^,    p-y^)      where 
-jS- — >  co.  We  may  fix  a  solution  by  fixing  the  height  at  co. 

In  Figure  3,  we  indicate  the  behavior  of  one  possible 
solution.  This  flow  has  been  computed  and  is  discussed  in 
Section  Ij.. 

We  muse  still  determine  \   from  (8)  which  we  rei^rite 

as 

du    _  dU_  +  cf dU 

>lds     '  Uds        IJ((1   _   u)2  m   2C2}    ds 

where 

u,c  * p — 7    ' 

s  =       r      /(&vr  +  (dcr    . 

TJ=1,C=0 
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2 
Along  the  curves  considered  C  /U  is  bounded  as  is 

-s ^ •  We  may  integrate  to  find  >x  for 

((1  -  U)^  -  2C^)ds 

Yl<   0 

log  (-  vt)  =  log  U  +  g(s)  +  log  Kx 

2 

where  g(0)  =  0  and  |f  =  2 — £IL and  K,   is 

as        U((l  -  U)    -   2CT)ds  x 

to  be   fixed.      Then 

(19)  f  \=  K-U  ee(s)   - 

At  U  =  1,  C  =  0  we  have  T\,=  k"   by  (l£)  and  (16)  and 


hence 


-k'1  =  K;L  . 


As  U  — - >  0  we  have 


H-  ±TT  ~  V    a        ° 


-u  ~  kx  e 

where   s  =  s(0,  0)  , 

For  \\  >   0  we  have  similarly 

log  \  =   log  (-  U)  +  g(s)  -  log  K2  , 
or 
(20)  Vt=   -K0U  e+S(s^   . 


Near  U  =  0,  Yy=  0  we  have 


-  ibC1-  k:1  e"S(s) 


. 
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r      -1 
If  u  =  r  U  =    U   is  to  be  continuous  on  >\j=   0  then  we 


must  have 


K2  e 


-S(s0)    T,-l  "S(so) 


=  Kx  e 


or 


K2  =  Kj_  =  -  k"1   . 


It  is  now  possible  to  determine  kp  .   The  line 


-1 


\   =  kl   corresponds  to  U  =  U, ,  C  =  C,   and  let  us  say 
Then 


s  =  s 


(21) 


1  ' 


"1       -1    SCSn) 

k2J"  =  k1J-  e     ux  • 


With  this  value  of  k0,  c?  is  determined  from  (18). 
To  summarize  lire  may  examine  the  solutions  in  the  (r,  t) 
plane,  see  Figure  2: 
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In  I  we  have  cavitation  u  =  0,  c  =  0.   In  II  the 
solution,  is  given  by  (3),  (19)  and  (20)  where  U  and  C 
are  given  by  one  of  the  solutions  in  the  U,  C-plane.   The 

state  III  is  again  at  rest  with  the  water  height  given  by 

L,2  ' 
6  2 

The  velocity  k„  of  the  reflected  shock  is  related  to 

k,   the  velocity  of  the  incident  shock  by  (21) 0 

3«   Finite  Depth  Problem. 

We  turn  now  to  the  problem  in  which  the  water  has  arbit- 
rary finite  depth  ahead  of  the  shock  but  the  shock  velocity 
becomes  infinite  at  the  origin.   Certain  conditions  will  be 
imposed  on  the  shock  path.  Let  the  path  of  the  incident 
shock  be  given  by 

(22)  r  =  R(t) 

or   t  =  H(r)   when   H(r)  <  0  . 

We  set 

dr 
(23) 

c   I  dii  C(k[FT 
t  d7 

in  the  region  k^  H(r)  >  t  >  H(r)   v/here  k2  is  given  by  (21) 
with  k^  =  1.   The  reflected  shock  is  given  by 
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(2k)  k2t  =  H(r)  , 

The  state  behind  the  reflected  shock  is  given  by 

(25)  u  =  0 


dR 
c  '  c2  dt 


First  let  us  show  that  (23)  is  an  sp  proximate  solution 
of  (2)  near  r  =  0  . 
With   \=  t/H(r) 


dH  -  TT 

dF  -  H  • 


(26)   ut+uur+2ccr  =  -y^[(-l-u)^Uv-   2CnC^-   U  +    (U2   +  2C2) 

(i  -  f )] 

_fl  +  Jfr  +  u     +  2  =  1[.  nn     +  2(l-U)C"1rlC   +   (3+^  *|t   ) 
c  c  r       r       tL      l  ■  \,  1  fir     §■ 

U  -   2]    , 

Comparing  (26)  with  (5)  we  see  that  the  differential 

equations  are  approximately  satisfied  if 
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(27)  |hh|  «  H2 

|H  -  Hr|  «  Hr 

These  conditions  are  satisfied  for  example  if  we  take 
as  shock  path  near  t  =  0, 

(28)  r  =  R(t)  =  t  log  (-  t) 
or  r  =  H(r)  log  (-  H(r))  • 


Then 


1  =  H(l  +  log  (-  H)) 


iL  -  i  =  H(l  -  H  log(-  H))  _ 
Hr         H  H  log(-  H) 


/.     log(-H)  \l+log(-H) 
\x   "  l+log(-H)j  logl-HT 


and 


7—^7 — rn-  — >  0  as  r  — >  0 
log(-  h) 


g-  fi2    i 

ri  "  ~  H   l+log(-  H) 


or  p  =  "  r+  log(-  H)  —  0   as   r  — >  0. 


At  the  shock  the  correct  conditions  would  be  given  by 
(12)  with  S  =  H"1  but  by  (23),  (22)  and  (16)  we  have 


(29)  uH  =  1 


• 


cH  =  0  ♦ 


Comparing  (29)  and  (llj.)  we  see  we  are  neglecting  terms 


■:.'. 
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•1/2 

of  order  H  '   .  Although  the  value  of  c   itself  at  the 

shock  is  wrong,  immediately  behind  the  shock  large  terms 

H  r 

of  the  order  —  which  behave  like  r-  enter. 

fit  t 

The  reflected  shock  is  given  by  t  =  kl  H(r)  .   In  the 
case  of  an  incident  shock  of  the  form  R  =  t   log  t   this 
means  that  the  reflected  shock  is  given  by 


r  =  k0t  log  k9t   k~t  log  t  . 


At  the  reflected  shock  x-xe  have  from  (23), 

k 
u  =  -^  U(kp) 

II 

k2 
c  =  -^  c(kP) 

H     * 

•-1 

and  shock  velocity  kpH   • 

The  state  immediately  behind  the  reflected  shock  is  then 
given  correctly  accoi'ding  to  (11),  in  (25): 

u  =  0 

c  =  c2/H  =  c2  ||  k"1  for  t  =  k^Hd?). 

In  the  whole  region  (III)  the  apparent  discrepancy  between 

(25)  and  (1)  can  be  justified  by  considering  the  neglected 

•  -1         • 
terms  as  being  of  order  R  where  R(t)  =  dR/dt.   For  ex- 
ample in  the  case  where  R(t)  =  t  log  t,   R   — >  0  as 
t  — >  0  . 
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In  fact  if  we  consider  solutions  of  the  form 
u  =  4{ij  U(R(t)/r) 


0  =  ^  0(R(t)/r) 


we  find  again  under  the  same  restrictions  on  the  shock  path 
that  U  and  C   satisfy  the  ordinary  differential  equations 
(5).   One  solution  corresponds  to  the  state  of  rest,  namely 
U  =  0.   Then 

4-4-'    R(t) 

C  =  constant.  — ~— - 
or         u  =  0 

c  =  constant.R 

as  in  (?.£)• 

It  can  also  be  shown  that  the  appropriate  perturbation 
problem  behind  the  reflected  shock  has  a  solution. 

h •      The  Computed  Flow . 

A  solution  of  (7)  issuing  from  (0,  1)  was  computed  as 
an  example  of  a  contracting  shock,,  Let  the  velocity  of  the 
contracting  shock  be  k,  .   The  solution  chosen  was  the  one 
for  which  the  height  at  infinity,  which  is  also  the  maximum 
height  of  the  incident  flow,  is  given  by 


Zmax  =klS   lim  „  Tt  ^(nj  =  1.59  g  *lc£. 


V. — >0 
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The  the  velocity  at  infinity  as  computed  is 

u   =  lira    ^l"1U(tt)  =  2912  k,  . 
00   >l— >  0  x 

The  velocity  of  the  reflected  shock  is  1„65  k   .   The  height 

-12     X   ' 
behind  the  reflected  shock  is   7»^8  g  k-,   or  q..7  times 

the  maximum  height  of  the  incident  flow. 

2    ? 
In  Figure  3,   gZ\  =  C   is  plotted  as  a  function  of 

ifW=  U  •   The  shock  polar  is  the  curve  given  by  (12). 

In  Figure  !+,  %^   is  given  as  a  function  of  U  , 


1 

■ 
... 
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solution 
of  (7) 


-1.20  -1.00  -.80  -.60  -,h0  -.20   0    ,20   .1+0   .60   .80   1.00 


Figure  3 
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1.1C11   -i.co-iao   -IbO   -ltj.0 
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-.60 
-.70 


.U.0    '  .60    '  .dO    '  1,00 


Figure   L\. 
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